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Using the frameworlt of the non-local light-cone expansion a systematic study is performed for 
the structure of the twist-2 contributions to the virtual Compton amplitude in polarized deep- 
inelastic non-forward scattering for general nucleon spin with an additional scalar meson in the final 
state. A useful kinematic parameterization allowing for appropriate triple- valued off-forward parton 
distribution amplitudes is given. One-variable amplitudes being adapted to the fixed parameters of 
the extended Bjorken region are introduced by decomposing the Compton amplitude into coUinear 
and non-coUinear components. These amplitudes obey Wandzura-Wilczek and Callan-Gross like 
relations. The evolution equations for all the distribution amplitudes are determined showing that 
the additional meson momentum does not appear in the evolution kernels. The generalization to n 
outgoing mesons is given. 
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I. INTRODUCTION 



Compton scattering of a virtual photon off a hadron, 

7*(<Zi) + H(pi)^7*((72)+H(p2), (1.1) 

is an important process in Quantum Chromodynamics. This general process covers a series of different reactions 
through which a variety of inclusive informations on the short-distance structure of nucleons become accessible at 
large space-like virtualities. It is also closely connected to the spin problem of the nucleon. The case of forward 
scattering pi = P2 = P describes deep inelastic scattering (DIS) off unpolarized or polarized targets which is widely 
discussed in the literature, see e.g. jl], ^, and pi ^ p2 corresponds to the generic (ordinary) non- forward virtual 
Compton scattering ||, ||, ^, 0, for recent reviews see In the special case g| ~ 0, i.e., when the outgoing photon 
is real, this process is called deeply virtual Compton scattering (DVCS). 

Experimental results on polarized and unpolarized (deeply virtual) Compton scattering were reported in ^ |l^, 
O, Il3|. The kinematic domain of some of these investigations is bound to rather low values of Q^. Experimentally 
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the final state in deep-inelastic non-forward scattering contains aside the (virtual) photon and final-state hadron, 
Eq. (1.1), a series of other hadrons, which even may emerge at the amplitude-level. The latter process is much more 
likely and of greater practical importance than that of a single isolated hadron -fl + Hi ^ 72 + H2, which was studied 
before I, I I ||. 

In this paper we extend the description given in the ordinary non-forward case in Refs. ^, Q to physical processes 
with an outgoing scalar meson, 



7*(gi) + H(pi 



7*(g2)+H(p2)+M(A;) 



(1.2) 



to investigate which of the properties derived in Refs. ^ remain valid in the more general situation and which 
are changing to account for more realistic experimental situations which allow for additional studies of distribution 
functions emerging in non-forward scattering [p7[. 



If one looks for a diagrammatic representation of the amplitude for the process (1^) then even in the special kine- 



matics of the Bjorken region for (p2 + k — pi) being small, there appear different production mechanisms. However, 



if one has in mind = i ^ 



and p2-k <^ M then soft processes between the two final particles and the 



incoming particle | pi ) are essential and we have to use a generalized distribution amplitude (p2, k\ 0\pi). For other 
cases one may try to make models which us e w ave functions of the nucleon ( p2 \ and the meson ( fc | . 
The Compton amplitude for the process (1^) is given by 



T^,{P+,P-,k;q)=i / d^xe- 



P2 



&;fc,0 



RT 



J, 



Pi,Si 



where 



P± = Pf±P,=P2±Pi + k, g = 1 (gi + q2) 



(1.3) 



(1.4) 



and k are chosen as independent kinematic variables. As usual pi (^2) and qi (52) denote the four-momenta of the 
incoming (outgoing) nucleons and photons, respectively, and S2 are the spins of these nuclcons and k is the 
momentum of the outgoing scalar meson. p± = p2 i pi denote the independent kinematic variables in the ordinary 
non-forward case. In principle, this set of kinematic variables can be extended to the case of n outgoing (scalar) 
mesons of momenta fc^, z = 1, . . . , 71 with Pf — P2 + X)i ^i- 

In ordinary non-forward Compton scattering the generalized Bjorken region is defined by the conditions 



keeping the variables 



V = qp^ 



and 



q 

qp^ 



and 



qp_ 

qp+ 



ql 



ql 



2v 



(1.5) 



(1.6) 



fixed. This definition has to be extended now by taking into account the additional momentum fc. We define the 
extended Bjorken region for light-cone dominated QCD processes by the conditions (^]^) keeping the following three 
variables 



i = - 



qP+ 



qP_ 

qP+ 



and 



X ■ 



qk 

qPl 



(1.7) 



fixed. This can be extended furthermore to the case of n mesons by the obvious generalization of P+ and keeping 
Xi = {qki)/{qP^) fixed. In the limit ki — > 0, these definitions reduce to the usual generalized Bjorken region. 

For completeness, we list the different kinematic domains for forward and general non-forward processes and the 
related scaling variables. These domains are distinguished as follows: 

• Bjorken region for forward scattering; fixed quantity: ^ (or xb = ^qi/^qi.pi) 

• Generalized Bjorken region for ordinary non-forward scattering, including DVCS kinematics for — 0; fixed 
quantities: ^ and 77. In the special case of DVCS S, — —r] holds. 

• Extended Bjorken region for non-forward scattering with a single outgoing scalar meson; fixed quantities: ^, rj 
and X- 

• n-Extended Bjorken region for non-forward scattering with n outgoing scalar mesons; fixed quantities: ^, rj and 
Xi, ■ ■ ■,Xn- 
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By conservation of momentum, 

=P2 + 92 + Ej^j> (1-8) 

one can show that ^ and 77 possess the same interpretation as in the ordinary non-forward case. Especially, also the 
relation ij = (^qf — gl) /2i' holds for processes including the outgoing mesons. The variables Xi describe the momentum 
fractions of the mesons ki in the infinite momentum frame defined by the momentum P_|_. In the following we r estri ct 



the consideration to the case of one additional final-state meson and return to the case of n- mesons in Section VII. 

It is important to remark that all physical processes mentioned above are distinguished only by taking different 
matrix elements of the same renormalized operator [^8| , namely the renormalized (i?) time-ordered (T) product 

f^.(x) ^ RT [j^ (I) J. (-|) S] , (1.9) 

where J^{x) —:'ip(x)j^'ip{x) : denotes the hadronic current and S is the renormalized S*— matrix. Near the light-cone, 
— > 0, this operator will be decomposed via the non-local operator product expansion into a series of non-local 
light-ray operators and the corresponding coefficient functions : 

Tp,v{x) J d\ Crf_cu{x^,Hix;^'^) RT[0^{kix,K2x)S] 

+ higher order terms . (1-10) 

The coefficient functions are singular on the light-cone. They are entire analytic functions with respect to KiX resulting 
in a restricted integration range — 1 < < 1. The unrenormalized light-ray operators in this expansion are given by 

0^{kix,K2x)= : '4> (kix) T U {kix, K2x) tp {k2x) : (l-H) 

with some specified F-structure of Dirac matrices and the usual path-ordered phase factor, 



U (kix, K2S:) 



P exp^-i5 J dri^A,, (t£)^ , (1.12) 



ensuring gauge invariance. Here is the gluon field, g the strong coupling constant and i is a light-like vector 
depending on x via a non-null subsidiary four-vector C, 



(xC) - V {xCf - x'^C ) with > 0. (1.13) 



In Eq. (I.IC) the flavor structure has been suppressed. Eventually, in the singlet case, also operators containing the 
gluon field strength F^i, and their dual F^^i, have to be taken into account. The contributions which contain four or 
more (anti)quark fields will be denoted by 'higher order terms', possibly also together with (powers of) the gluon field 
strength, etc. By construction, the non-local light-cone expansion, depending on the order of terms being taken into 
account, leads to a (sub)asymptotically relevant part and a well-defined remainder being less singular, see Ref. 

Taking matrix elements of the operators and performing a Fourier transformation in the expressi on ( |1.3| ) leads 
to the physically interesting non-perturbative distribution amplitudes. Because the operational input ( |l.9| ) for these 
amplitudes is the same for the different processes described above, the evolution equations of these amplitudes are 
determined by the renormalization group equations of the operators and, therefore, also by the same anomalous 
dimensions. 



This paper is organized as follows. In Section^ we discuss the quark-antiquark operator (2.13) as operational input 



for the (extended) Compton amplitude. This allows one to use the known anomalous dimensions of this operator to 
write down the evolution equations for the relevant distribution amplitudes used in the process including an outgoing 
meson. The evolution kernel required for these evolution equations is determined by the anomalous dimensions of 
the quark-antiquark operator and further computations of Feynman diagrams are not necessary, at least at one-loop 
order. Furthermore this operator possesses a known twist decomposition given in Refs. ||l^, |l8|. In this paper we only 
determine the twist-2 part of the Compton amplitude in the simply extended Bjorken region which generalizes the 
twist-2 representation known for forward and ordinary non-forward scattering. 



In Section [II the twist decomposition of the matrix elements and the necessary decomposition into suitable kine- 
matic factors is performed thereby using the hadron equations of motion. This section also includes the definition of 
the distribution amplitudes used in the process considered. 
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In Section IV the twist-2 part of the Compton amplitude is calculated. It will be shown that, in the extended 



kinematic region, the amplitude depends on the three scaling variables ^, rj and x ^^^d is given by triple-valued 
distributions. 

In Section ^ we extract integral relations contained in the Compton amplitude and Section VI includes the evolution 
equations obeyed by the distribution amplitudes. 

In Section VII we add some remarks on generic properties of the results obtained in the preceding sections and 
consider the case of n outgoing mesons; Section VIII contains the conclusions. 

Appendix A contains the projections of the twist~2 operators on the light cone. In the Appendix ^ we construct 
the helicity basis for both photons and calculate all the helicity amplitudes. These projections show in explicit form 
that the process (|l.2|) is current conserving on the level of the ^-matrix. 



II. OPERATOR STRUCTURE 



In this section we discuss the operator structure of the Compton amplitude which will be used in the following sections. 
For brevity we discuss only the flavor non-singlet case and drop all flavor structures in the operators. The construction 
for the flavor singlet case is to be carried out similarly. 

As was mentioned above, the operator input for the Compton amplitude T^^[P±^ k\ q) is given by the renormalized 



time-ordered product of two electromagnetic currents, Eq. (1.9). Obviously, in lowest order in the coupling constant the 
S'-matrix can be set equal to one. Then, applying the Wick-Theorem to this time-ordered product and approximating 
the quark propagator near the light-cone by 



S{x,m) — (i ^ + m) A (x, m) 



(2.1) 



m 



2^2 (^2 _ ig)2 (3.2 _ 

one obtains the following description of the operator T^i,{x) 

1 



— (2i + m^) + 



1287r- 



■ln(-x^ -fie) (Si + m:^) 



T 



■'»l|)^"(-f 



1287r2 



■ In i-x^ 



27r2 (a;2 - ie)^ 87r2 {x^ - ie) 

X (: (a;/2) l^^lali- ^ i-x/2) : - -.i) (-.t/2) 7^707^ V {x/2) 



(2.2) 



IGtt' 



■ In (-a;2 + ie) 



47r2 (.t2 - ie) 

(: {x/2) (7^^ V {~x/2) ■.-■.4, {~x/2) V {x/2) :) 



+ : {x/2) ^ {-x/2) : + : {-x/2) ^ {x/2) 



+ ■.^{x/2)-f^^p{x/2)^p{-x/2)-f,^jj{-x/2) : + 



^4 (x2 - ie)" 



The renormalization symbol R is suppressed in the notation. The four-quark operator ip^fj^tpip^^ip has no singular 
coefficient on the light-cone and is of minimal twist-4 and the last term stems from disconnected diagrams. Both 
terms are therefore discarded, see Ref. Q. 

The quark mass terms resulting from the mass dependent part of the operator (i^ 4- m) are less singular and will 
also be omitted. In fact, the lowest twist contribution of all quark mass terms is contained in the operator 



i(:i^(x/2) iaf,,t^{-x/2): - ■.^{-x/2) ia^, V (2;/2) :) 



(2.3) 



with (7^^ = (i/2)[7^, 7^] and is of t wist -3. 

The first term of the expansion (2.2) is of main importance, because it contains the leading light-cone singularity 
and its minimal twist contribution is of twist-2. It also contains terms of higher twist (trace terms) and quark mass 
terms resulting from the mass dependence of the scalar propagator A {x.m). These terms are also less singular on the 
light-cone. 
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Using the standard relations 



(2.4) 
(2.5) 



where we indicate that S^avp is symmetric in and a/?, and noticing that the operator (2.2) is dominated by the 
Hght-cone singularity, w 0, one arrives at 



2 I '-'fJ^au/i 



X X 

2'~2 



S 



O 



5/3 



X X 

2'~2 



TT^ {x"^ — ie) 

with the (unrenormalized) light-cone operators 

(i/2, -i/2) = ^ (: ^ (5/2) 7^ ^ (-i/2) : - : ^ (-5/2) ^ (i/2) :) , 
0^''(i/2,-i/2) = i(: i/i(i/2) 7VV'(-5/2) : + : V^(-5/2) 757'^V(£/2) :) 
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(2.6) 



(2.7) 
(2.8) 



This expansion has to be viewed as a simple form of the non-local light-cone expansion (1.10) from which the suitable 
F-structures can be read off. Considering general gauges also the phase factors U (kiX, K2x) must be included into the 
light-cone operators and O^^ . Taking into account also higher order terms of the S'-matrix additional operator 
structures come into the play which, of course, are sub-leading. 

The K-integration appearing in the light-cone expansion (1.10) connects the coefficient functions, given here at Born 
level, 



Ci/2 (a;^,Ki,K2) = 



27r2 (x2 - ie)^ 



(2.9) 



with the respective operators. As is known from the general analysis of the light-cone expansion the Fourier 
transforms of the coefficient functions in Eq. ( 1.10| ) are entire analytic functions with respect to the variables xpi 
which leads to a restriction of the integration range of ki and K2 to the interval [—1, +1]. 
For later convenience we introduce the variables 



K± — ^ (k2 ± Hi) with Kl,2 = K4. T 

and use the following integration measure in the light-cone expansion 
D^K = dKi dK2 6* (1 - Ki) 6* (1 4- Ki) 6* (1 - K2) 6* (1 + ^2) 

= 2dK+ dK_ e{l- K++ K-) e{l + K+- K-) 9{1~ K+- K-) 6 {1 + K+ + K-) 

In terms of the variables k+ and k_ the coefficient functions read: 

i 



(2.10) 



Cl/2 = 



6 (k+) 5{k- + ^)tS (k- - i) 



47r2 (x2 - ieY 

Using these conventions one obtains the following representation of the operator T^^{x): 



RT 



J, 



(f)M-f)^ 



Ci (a;^, K+, K_) S^j,jj\°'^ Xa Of3{{K+ - K_) (k+ + k_) x) 

- 16*2 {x^, K+,K^) 6^^°'^ Xa 0^{k+ - K_)x, (k+ + K^) x) 



(2.11) 



(2.12) 



In general, the renormalized non-local operators Op {kix, K2x) and (kiz, K2x) containing the phase factors 
U (kix, K2x) are given by 

Ofj {kix, K2x) = ^RT (: 4' ('*ii) 7/3 U {nix, Kax) ■(/; {hi2x) : — : -0 {k2x) 7^ U [nix, K2x) ip (kix) :) S (2-13) 

0^(Kii,K2x) = ^RT (: 7/) (kix) 7^7^ t/ (ki£, K2i) ('«2i) : + : ■0 (^2^) 7^7/3 ('tis, ^2^) V' (kiS) • (2-14) 



In the following we assume that the operators O/3 and Og are renormalized quantities. 
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III. TWIST DECOMPOSITION AND MATRIX ELEMENTS 



The non-local operators Op {kix, K2x) and O^^ {kix, K2x), Eqs. ( ^.IS ) and (2.14), contain contributions of different 
twist. Here, the notion of twist is used in its original form [M as 



(geometric) twist (r) = scale dimension (d) — Lorentz spin (j) 



(3.1) 



The operators appearing in the expansion (2.12) of the Conipton amplitude have to be decomposed into their various 
twist parts. On the light-cone they contain contributions of twist-2, 3 and 4 



/-j5 tw2 



/-jtw3 



/-jtw4 



/-)5 tw4 



(3.2) 
(3.3) 



however, off the light-cone they contain an infinite series of growing twist. 

A group theoretical procedure of the twist decomposition has been worked out in Refs. |p7|, |2C | and applied to 
various physically relevant light-ray operators. As a result, the twist-2 part of the operators (2.13) and (2.14) can 
be constructed out of the twist-2 (pseudo) scalar operators o(5)tw2 (^i^^x, K2x) = x^O^^ (nix, K2x) by applying the 
interior derivative 



(l + xd^ dp - ^xp 8P 



with 



d^ 



, 



(3.4) 



on the un-decomposed light-cone operators and performing a subsequent r— integration (which stems from the nor- 
malization of the local operators): 



O^^"^ {kiX , K2X) = dr luT dp O {kiTX, K2TX) 

Jq 

O {kiTX, K2Tx) 



(3.5) 



dr 



dp 



-hiT xp 



According to its structure the tracelessness of the operator ( p.5[ ) which corresponds on the light-cone to the requirement 
= 0, is trivially fulfilled due to the property of d'^. An analogous relation holds for the axial vector 
and pseudo scalar operator. The twist-3 and twist-4 parts, however, cannot be constructed out of the (pseudo) scalar 

operator since the latter, when restricted to the light-cone, is already of twist-2. 

Since we want to extract the twist-2 part of the Compton amplitude, relation (3.5) will be applied to the matrix 
elements of the operators considered : 



{p2, k\ O^p"^ (kix, K2x) \pi) 



dr 



dp + - Inrxp d"^ 



(p2: k\ O {kiTX, K2Tx) \pi). 



(3.6) 



Here, the spins of the nucleons have been suppressed in the notation. Let us mention that the geometric twist 
decomposition of the matrix elements is due to the twist decomposition of the non-local operators. Usually, in 
phenomenological considerations another notion of twist, called 'dynamical' twist, is considered which has been 
introduced in the decomposition of the (forward) matrix elements by Jaffe and Ji ||2l|] . The interrelation of geometric 
and dynamic twist was considered in Ref . ||2^ . In the case of lowest twist-2 there appears no difference, but for higher 
twist the mismatch of dynamical twist with respect to geometric twist leads to differing structures. 
Because of translation invariance, 



{p2,k\ 0{kiX,K2x) \pi) 



AK-\- xP- 



(p2, k\ O {~K^X, K^x) \pi) , 



(3.7) 



it is more convenient to discuss the centered operator O {—kx, kx). Henceforth, for brevity, k_ will be denoted by k. 

In the ordinary non-forward case one usually parameterizes the scalar matrix element by a Dirac- and a Pauli-type 
contribution, 



/Ci(5,p2,Pi) = 5^ w(p2) 7/3u(pi) , 



1 



/C2(i,P2,Pl) = x'^ U{p2) <JpaP°Lu{pi) , 

TOO 



(3.8) 
(3.9) 
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where u{pi) and u{p2) are on-shell spinors of the incoming and outgoing nucleons, toq is a dimensional mass scale 
which is kept fixed, and a^a — (i/2) [7/3 , 7a]- Because of p_ ^ the Pauli type factor IC2 vanishes in the forward 
limit. Using these kinematic factors the matrix element can be parameterized as follows: 

2 

{p2\0{-KX,Kx)\pi) = y^/Cg {x,P2,Pl) fa {kXP+, KXP-, PiPj, fM^) . (3.10) 

a=l 

Here, the coefficient functions fa are the distribution amplitudes in a;— space. They depend on Kxpi and all possible 
products of the two momenta, pipj — {pi,P27PiP2} , as well as on the renormalization scale /i and the coupling 
constant g. 

If an outgoing meson is present in the process, which means that one has to contruct the matrix elements 
^^25^1 0{—Kx,Kx) the representation ( |3.10 ) of the scalar matrix element must be modified because of the 
presence of the additional momentum k. Especially, further kinematic structures occur. They can be determined in 
a straightforward way by using the following parameterization of the scalar matrix element: 

{P2,k\ O {-KX,KX) \pi) = U {JP2) A{KX,P2,k,Pi) U (pi) , (3.11) 

with 

A («i,p2, k,pi)^AI + B„ 7" + C[^i3] t"" + Da 7V ■ (3.12) 

It is an easy task to find the general structure of A — D allowed by the momenta pia , P2a , ka , the metric gap and the 
Levi-Civita tensor epa-yS, demanding A not to be a pseudo scalar. 

Using the equations of motion for the hadronic momenta pi^ and p2f_i with M denoting the nucleon mass, 

7Va "(Pi) = Mu{pi) , (3.13) 

U{P2) 7"P2a = Muip2) , (3.14) 



the decomposition of the matrix element (3.11) is : 

IC[ = {uu) , (3.15) 
JC2 = {u$ u) , 

/C3 = — (-U ^ u) , 
mo 

/C4 = (uiji l/l: u) , 

mo 

/Cg = iepa'jS x'^ P2 k'^ pI {u-y^ u) , 
m^ 

where an auxiliary mass mo has been introduced in order to get kinematic structures of equal dimensionality. Then, 
the scalar matrix element is parameterized by a sum over the five kinematic factors /C^ as follows: 

5 

{p2,k\0{-Kx,Kx)\pi) = ^/Cl(£,p) /g (rap, PiPj, /i^) , (3.16) 

where p = {pi} ~ {pi,p2, k} generically denotes the multi-vect or in the space of all the (three) hadronic momenta. 

Although one possible set of kinematic factors is given by ( [3.16| ), it will be more convenient to choose another 
one which is also linearly independent and contains the original kinematic factors /C^. Mimicking the Dirac- and 
Pauli-structures we choose 

K.i{x,p2,k,pi) = x'^{ujpu), (3.17) 



1 



u 



/C2 (i,P2,fc,Pl) = 5 [uCJpaP- 

mo 

/C3 (i,p2,A:,pi) = x'^ —{ukpjak"' u) , 



m, 



K.4{x,p2,k,pi) = x'^ —{uapak"' u) , 
niQ 

JC^{x,p2jk,pi) = x^^ —ispa-iSP'^k"' Piiui^ u) , 
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for the matrix element of the scalar operator and 



JCl{x,p2,k,pi) = i^iuj^-fpu) , (3.18) 
mo 

mo 

^l{x,P2,k,pi) = x'^ — gie/sa^aPa u) , 



for the matrix element of the pseudo scalar operator. For explicit calculations it is important to note that each factor 
id^a^ has a linear x-dependence, 



K.a{,x,Tp) = xf^ ICapip) , 
Using the equations of motion again, the factors )Ca equal the following combinations of /C^ 





= 1^2 ' 






IC2 


= -( 




mo 






IC3 






mo 


IC4 


= ^( 




mo 


IC5 





(3.19) 



M 

'r xk + xpi) IC'i — 2i 

mo 

-1 

-3 ' 



This shows that the /Cq constitute a suitable set of kinematic factors for the scalar matrix element. In the limit fc ^ 
they obviously reduce to the Dirac- and Pauli-structures. Analogous statements hold for Ki\. 
The decomposition of the matrix element of O (— kx, kx) and {—kx, kx) now reads 

5 

{p2,k\ O {-KX,Kx)\pi) = ^JCa{x,p) fa{KXp, PiPj, fi^) , (3.20) 

a=l 
5 

{P2:k\ {-Ki,Ki)\pi) = ^/C^ (i,p) (Kip, PjPj, . (3.21) 



In principle there is also a dependence on variables like PiPj x^. Because the latter dependence vanishes on the light- 
cone, it will not be discussed in the further considerations. As has been shown in Ref. |l^ the whole x^— dependence 
is governed by harmonic extension off the light-cone if the operator structure is already given on-cone. The 
/x^— dependence is governed by the evolution equations obeyed by the functions fa and will be discussed in Sec- 



tion VI 



As next step, a Fourier transformation of the functions fa is performed, 

/„(Kip)= / D^z e-'^^tP^'Valz) . (3.22) 



Because also the functions fa (ftxp) are entire analytic in the variables Kxp, the support of their Fourier transforms 
fa (z) is restricted to the interval [—1, +1] in the variables z = (zi, Z2, z^). Therefore, the measure 

D^z = dzi dz2 dz3 9{1- zi) 9{1 + zi) 6(1- za) 6* (1 + za) 6* (1 - Z3) 9 {1 + zg) (3.23) 



has been introduced to realize this support, pz is simply the product of the vectors p and z, see (3.24). To get a 
representation in the momenta P± and k one introduces the variables 

1 / ^ 1 / 

Z+ = -{Zl+Z2j, Z_ = -(Z2-Zl), Zfc = Z3-Z2, 
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with 



Using the abbreviation 



P+z++P^z^ +kzk 

Pi Zl+P2Z2 + kz3 = pz , 



(3.24) 



the scalar matrix element is represented by 

(p2i fc| O {~KX, KX 



)\pi) = ^/^^(i,?) / D^z e--^^/,( 

a=l 



Z+,Z-,Zk) 



(3.25) 



In the following th e exp licit summation over a wi ll be omitted, but will be indicated by the position of the index a. 
The expression ( [3.25| ) will be inserted into (3^) to calculate the matrix element of Oj^^. Thereby, it is important 



and Zk Zk/r we get the following form of the matrix element, for the general case ki ^ —K2- 



to state that the r— scaling in equation (3.6) refers to n and not to a;. Performing a change of variables z± z±/t 

(3.26) 



(p2, k I O^f^'^ {kix, K2x) I Pi) 



dz^ dz_ dzk 



dr 



9/3 + - In r Xf3d 



X e 



\ T T T / 



where the abbreviation 



8 (r, z+, z_, Zfc) = 2 6* (r — z+ + z_) 6* (r + z+ — z_) 6* (r — z+ — z_) (t + z+ + z_) 
x9 {t — z+ — z^ — Zk) 6 {t + z+ + z^ + Zk) 

has been used. Carrying out the differentiations we get 

^1 dr 



(3.27) 



{p2,k\0'^^K,£,K2i)\pi)= f dz+dz_dZkf ^e-('^+-^---^)/C"''(p)/a(^,— ,-) 9 (t, Z+ , Z_ , Zfc) 

Jr;' Jo t \ t t t / 



gfSp + i [h+tP-p - nVp) Xp - — Inr ( {k+tP^ - nVY Xp - 2i {n+rP^p 



n2 ~ 



(3.28) 



This form of the matrix element is yet rather complicated. Since only the centered operator is needed in the following 
considerations, we set k+ = 0. 

To derive a simpler representation for the matrix element, the r- integration will now be comprised into the functions 
Fa and ^ the latter denoting the trace part, defined by 



1 



Fa(z) EE (z+,Z_,Zfe) = / dr — /a , ,— , 



Z+ Z_ Zk 



In^ 



T T T 



Z+ 2- Zk 



Fr(z) EE i^-(z+,z_,zfc)= dr^/,^,— ,^ , 



T r r 



with fa given by 



f fZ+ Z^ Zk\ _ fZ+ Z_ Zk\ , 

Ja[ , , — ]=Ja[ , ,— B T, 

\ T T T / \ T T T / 



Zk) 



(3.29) 
(3.30) 

(3.31) 



Obviously, these distribution amplitudes are not independent, and the restricted integration range in z— space finally 
is contained in the support properties of t he dis tribution amplitudes Fa (z^, z-, Zk) and i^*"^ (z^, z-, Zk)- 
After the substitution of Fa and Fa in ( 3.28| ) with k+ = one obtains 



(P2, k I Oj}''^ {-KX, KX) I pi) 

dz+dz_dzfc e^'''^^ 



(3.32) 



{gpp - iKVf3ip)Fa{2)~^- {n^V^ip + 2iKVp) F'Jiz) 



IC^^ip) 
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for the matrix element of the twist-2 part of the vector operator O^. The matrix element of the axial vector operator 
possesses a similar representation 



tw2 



dzj^ dz- dzk e 



-IK XV 



Xf3 



(3.33) 



/C5°''(p) 



with and i^^'"' defined analogously to Fa and Fl\ by exchanging fa and fl in ( ^.29| ) and ( ^.30| ). 

Here, some general remarks are in order. First, the triple-valued distribution amplitudes i^i^'' (z+, z_, Zfe) and 

^a^'*' (•2+,2-,Zfc) are uniquely related to the twist-2 (axial) vector operators and, in principle, should have been 
marked by the related twist t — 2. However, since we do not consider higher twist this has been omitted. Second, 
every distribution amplitude of definite twist - also off the light-cone - depends only on the z— variables and, possibly, 
on the momenta p. The x— dependence is completely contained in the accompanying factors including, of course, 
the exponential e""^^'''. These general expressions which are given in terms of Bessel functions of the argument 
{k / 2) {xVy — x^V-^ have been determined in Ref. |l8j (For the twist-2 ca se of PIS t his st atement has already been 
made in Ref. [^.) Restricting onto the light-cone leads to the expressions ( 3.32| ) and (3.33). Third, these properties 
hold for the operators of definite twist and are transposed to the corresponding matrix elements, independently 
how many particles (momenta) occur in the incoming and outgoing states. Therefore, if a Fourier transformation 
containing these matrix elements has to be performed this can be done by simply replacing x ^ x. (In Appendix A 
we give these expressions explicitly to geth er wit h the ir restriction onto the light-cone.) This will be applied in the 
next Section for the expressions Eqs. ([3.32) and (3.33). 



IV. COMPTON AMPLITUDE 



Now, we are in a position to compute the twist-2 part of the Compton amplitude (1.3). Of course, we need it in the 
extended Bjorken region and, therefore, can restrict our consideration to the neighborhood of the light-cone. Thus, 
we will take the matrix elements ( [3.32 ) and (3.33) with x replaced by x. Merging every thing toget her, we use the 
representation ([2.12|) for the time-ordered product of two hadronic currents, insert (3.32) and (3.33) for the matrix 



elements of Of'''' and O^^""^ and obtain: 

T^Z^{p2,k,p,,q) = i J e"^- X^ Cl{x\K+,K)Sp,^''^(^P2,k\0f^{{K+-K)x,{K+ + K)x)\p,) (4.1) 

-iC2{x^,K+,K)ep^°''^(^P2,k OIi''"'^{{k+^ k)x,{k+ + k)x) pij) . 

The arguments and k of the coefficient functions Ci (x^, k+, k) are fixed at k+ = and k = ±1/2, so that we can 
use the symmetry properties of the operators Ofj and 0|, 



Ofj{ — KX,Kx) = —Of){KX,—Kx) , 

0^{—Kx,Kx) — 0^{kx,—kx) . 



(4.2) 



For the computation of the Compton amplitude it suffices to know the operators at these given points, but for the 
investigation of the evolution of the matrix elements their representation at general values of and k is needed, see 
Section VI for the details. 



Performing the ^-integration one obtains 



TlZ^{p2,k,pi,q) 



dz+dz_dz, /— e ^ 



iQx 



(4.3) 



a(3 



90P - \Vfi xp\ Fa - xp (^V^ xp + '-Vp ) Fl' 



VbxA F^JC'-p 



where the abbreviation 



Q = q-^V^q-^ {P+Z+ + P_z^+k Zk) 



(4.4) 
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has been used. The i^-term which results from the traces of the twist-2 operator vanishes for the axial matrix element 
because XaXp is symmetric. As a last step in the computation of the Compton amplitude the Fourier transformation 
is carried out by using 



d^ 












d^ 












d^ 






2^ 




{x^ — lef' 



+ le 

9al3 



2i- 



QaQ/3 



(Q2 



= 2 



Sq/J Qp + 5/3p 2a + 5ap S/3 „ Qa 2/3 2p 



(22 + ier 



and the summation over a and /3 is performed in the symmetric part of T"*^ by using the form (|2^) for 5^^"'^. Then 



for the Compton amplitude we get the result 

1 



T. 



tw2 



dz_|_ dz_ dzj. 



22 



-1£ 



a/3 „ 



5/3p 



^/3 2p 



5pp (gz. - ^z.) + 9vp (gp - T'p) - g^,y [qp - Vp) + 
1 



22 



^5 ^5ap 

QpVy + Q,V^,-gp,VQ 



(4.5) 



dpf'Pp 



QpV 



(2 2p2. - 5p.2')7'p - i(g.p2p + gpp2. - 25^,2^)7'^ 
^(2 2p2.-5p.22) I i^f/C^" 



which is expanded with respect to the functions , Fa and i^*'' with the trace term separated from the antisymmetric 
and remaining symmetric part. 

This structure of the Compton amplitude is a generic one because it is also valid for the ordinary non-forward case: 
The additional structures arising due to the momentum k are hidden in the summation over the structure functions 
Fa, Fa and Fa and in the definitions of V and Q. The reason for that result is an outcome of the twist structure of 
the operator which is the same for all the matrix elements under consideration. It holds also for the case of n outgoing 
mesons. 

For V = p+zj^ +P-Z- and summation over the Dirac- and Pauli-structures one reproduces the form of the Compton 
amplitude given in Ref. |Q. Here, the additional terms containing the functions F^'' arise, because the trace term in 
(3.6) has been taken into account. If k is not present in the process this term vanishes in the limit 



l^Pip u{pi) 
u{p 



2) Y'P2p 



(4.6) 
(4.7) 
(4.8) 



setting M and t = (p2 — Pi)^ to zero. Therefore these terms have been omitted in Ref. 0. If the five kinematic factors 
lCa{x,p) containing the momentum k are present this is a priori no longer the case since non-vanishing contractions 
like {u^t u) appear. Therefore, this term has been taken into full account here. 

In general, terms containing the product VplC"''' do not vanish in the massless limit, while terms containing V'^, 
PiPj, are small compared to the large invariants. This leads to the approximation 

(4.9) 



22 



— q.V 



qP+ C + (z+ + rjz- + xzk) - ie 



We apply these approximations to the Compton amplitude and get a representation for its twist-2 part in the massless 
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limit and in the extended Bjorken region, which depends exphcitly on the three variables z± and z^: 

1 1 



I dz+dz-dzk 



qP+ ^ + (z+ + 77z_ + xzfc) - ie 

ffMP (9"^ - 'Pi^) + 9vp [qp, - Vf,) - gf,,, [qp - Vp) ■ Fa (z+,z^,Zk) K."'' 



(4.10) 



+ g^^^Vp ■ F'J {z+,z-,zk) JC^'f - ie^."^ go gpp ■ (z+, z.,zk) K^"' 
1 1 



{qP+f (e + (z+ + + x^fc) - ie)' 



g^^- + - 'P^'P'^ - 9p.u q-V\ Qp ■ Fa {z+,z^,Zk) /C"" 
f (2 QpQ. ~ gp, q^ + gp. q.V) Vp ■ F'J {z+,z-,zk) K.'^" 

-ie^.^NaT'^jQp ■Fliz+.z-.zk) K.^'^P^ 



This form of the Compton amplitude will be used in the further considerations. 



V. INTEGRAL RELATIONS 



The variables z± and Zk are not directly measurable because they appear as Fourier variables of the distribution 
amplitudes fa- The scaling variable ^, however, can be regarded as a physical quantity and it is therefore quite 
natural to use the new variable 



t = z+ + rjz^ + xzk 



(5.1) 



as integration varia ble in the denominators of Eq. ( 4.10| ). 

By the definition (3.24) the vector V contains the variables z± and Z}~ and must be rewritten using the set {t, Z-, Zk}- 
This is simply done by introducing the combinations 



which leads to the representation 



TT = P--f]P+ , 
TT = k- XP+ : 

V = P+t + TTZ-+TTZk 



(5.2) 
(5.3) 

(5.4) 



The 4-vectors tt and tt define off-collinear directions w.r.t. to the direction P+; vectors along this momentum are 
denoted as collinear. Note that these vectors are non-forward still, since rj ^ 0. V contains collinear and off-collinear 
contributions, the former of which are associated with the scaling variable t only. It will turn out that these collinear 
parts play the dominant role in the process considered. 

Different powers of V contained in (4.10), of course, will lead to a whole series of structure functions corresponding 
to different moments in z_ and Zk- Therefore, let us generally define double moments of the triple- valued distribution 
amplitudes leading to single-valued distributions as follows; 



Fin, it;v,x) 



1 



1 dr 



dz_ / dzk z"' Zfe' F" {t - ?7Z_ - xZk, z^, Zk) 

t 







ni+n2 fa 



^ ' f nin2 \ X 



sign(t) 

— A-"^-"- • /„^„, (A; ry, x) 



with 



fin, (z-,V,x] ^ [ dz- / dzk zj r ^ * 



rjz- - XZk, Zk 



(5.5) 



(5.6) 



(5.7) 
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In fact, the values n.^ = 0, 1 occur in the antisymmetric part and the values rii = 0, 1, 2 in the symmetric part. To 
keep the discussion short, we consider the antisymmetric part of the Compton amplitude in full detail and give the 
result for the symmetric part only for the leading terms, i.e., suppressing the trace terms. The explicit calculation 
shows that the latter terms do not contribute to the leading order in v. 

Because the contractions VpJC"''' are only prese nt for the Dirac structure, K}'' = {u^''u), we discuss this term 
separately. Then, the antisymmetric part of ( 4.1C ) is given by 



T-itw2 



2i £u^ 



2i Saiy 



a/3 

qP+ 

a(3 la 



dz_|_ dz_ dzfe 



9Pp 



1p 



dt 



+ J-i 



i + t-ie qP+(^ + t- ie) 
9Pp F,%^ it) q, t {P+pPoS it) + ^pF!o it) + ^pFp'i it) 



i^ + t-ie) 



F!{z+,z.,Zk)IC"'''{p) 



(5.8) 



At this point it is necessary to perform a partial integration in the second term proportional to Qp using the formula 



dt 

^1 (e + < - ie) 
which leads to 



it;ri,x) 



dt 

-1 e + t-ie 



m-Fn^n, it;V,x)-t- 



fnin2 iti '?! X) 



(5.9) 



ptw2 



9- aP qa 

qP+ 



9l3pFi 



7 - 

-1 i + t~ie 



5a 
00 



qP+ 



p5a 
^00 



f5a 

loo 



5a 
10 



J 10 
t 



+ *rAFl 



'5a 



f 5a 
JOl 

t 



(5.10) 



^bap 



Several tensor structures contribute. Note that in the foregoing discussion no assumption has been made on the 
direction of the nucleon spin. As the polarization of the initial state nucleons in experiment is performed in outer 
magnetic fields, the direction of the nucleon spin is not related to other vectors in the system except for the condition 
Si.pi = to hold. 

The form (5.10) of the polarized Compton amplitude is very interesting because it includes a Wandzura-Wilczek 
(WW) like relation between the distribution amplitudes being associated to two of the tensor structures |3^. This 
relation becomes obvious, once the definitions 



G? it;v.x) 
(i;r/,x) 

Gg it;v.x) 

G^ (t;ry,x) 



^ floit;v.x) 

^ -f^Sit;ri,x)^ 

^ f!Sit;v,x) 
t 

^ f^Ht;v,x) 
t 



«ign(t) ^ 

T 

sign(t) ^ 
A2 



/oo (A;'7,X) 



f?Si^;v,x) 



^/o\"(A;^,x) 



(5.11) 
(5.12) 

(5.13) 

(5.14) 



are made. This leads to a very simple form of the antisymmetric part of the Compton amplitude, namely. 



rTitw2 _ r,- a/3 qa 



d. ' 

-1 e + t-ie 



5/3p(G? + G^) 



qP+ 



j^5ap 



All the above functions GJJ are expectation values of twist-2 operators. By definition G" and G2 obey the following 
integral relation 



G^ (t;r;,x) = -G? (i;ry,x) + ^ 



sign(t) ^ 

T 



G? iX;v,x) 



(5.15) 



This relation between two twist-2 quantities can be viewed as a generalization of the WW-relation known from forward 
scattering , . It is a property of the coUinear part of the polarized Compton amplitude Tj'^'j and connected to 
the 00- moment in z- and Zk- The off-collinear vectors tt and tt are connected to the new parton distributions G3 and 
G4, which obey an integral representation in terms of the functions /fp and /gf respectively, similar to G2. However, 
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unlike the case for Gi the respective functions do not appear in the Compton ampUtude. Therefore we obtain at the 
present level only one Wandzura-Wilczek like relation between the twist-two parts of the respective amplitudes. Of 
course all the functions Gi do receive higher twist contributions, which were not discussed in the present paper. Also 
these contributions as emerging in the different amplitudes may obey similar integral relatio ns. The generalization 
(5.15) of the WW-relation has been obtained in Ref. |^ for the ordinary non- forward process (1.1) without outgoing 
meson. The foregoing discussion shows that it remains valid for the more general process (1.2). 



The Wandzura-Wilczek relation Eq. (5.15) is an identity between physical amplitudes which emerges at the level of 
geometric twist-2 as a consequence of the fact that the distribution f^^it, rj; x) appears two times in the decomposition 
of the Compton amplitude. These relations which determine the (geometric) twist-2 content of the dynamical twist-3 
distributions have to be called geometric WW relations |2^, ^ ^ . Th ey a re obtained by using solely group 

theoretical means lying behind the definition of (genuine) geometric twist, Eq. ( ^.l[ ). These WW-relations have to be 
distinguished from the so-called dynamical WW-relations being obtained by using the QCD equations of motion as 
has been done by |29). To bring it to the point: Geometric WW relations are written for distributions having equal 
geometric twist, whereas dynamic WW relations are written between distributions of equal dynamic twist. Despite 
having the same formal structure their physical content is different. For a detailed discussion of these relations in 
the case of meson wave functions see Ref. |27]| which, however, with appropriate modifications also holds for general 
non-forward amplitudes (See also Ref. |2^ where the case of quark distributions is discussed). 

We now return to the discussion of the term containing the contractions VpJC""^. To begin with, we first remind 
that only the Dirac structure K.^^ is of relevance in this case and one may use the approximation 



t+{l-^) z_ + (1 - x) ^fel /C(5) 1" 



(5.16) 



in the limit of vanishing nucleon masses. Even more general, one can state that the product qpK^^^ is of order v 
and VpKP'^ is of order /i with /i ^ z^. The additional terms are therefore of non-leading type, but are interesting 
because these structures arise due to the meson momentum k. 

Inserting the above approximation and performing again partial integrations leads to the complete form of T'j^^. 



2ie "'^ — 



dt 



£, + t-ie 



G2) 



qP+ 



{P+13 G2 + 7r/3 G3 + 7f/3 G4) 



IC^^P (5.17) 



qP+ 



{P+P G2 + TT^ G3 -t- TTf3 



5 Ip 



with 



g; 



g; = t 



Gi 



-^00 



f51 

Joo 



/•51 
J 10 



f51 
JlO 

2t 



2^2 



+ i^~x){F, 



51 



^ +a-v) F^o'-^ +il-x) F^t-^ 



2t 



+ (i-x) 



p51 
^02 




f51 
J 02 



(5.18) 
(5.19) 
(5.20) 



The meson momentum k and the momentum transfer q are connected to similar structures, but the functions G^ 
related to k are far more complicated than G^. However, despite that fact the expressions within the parentheses in 
Eqs. ( 5.1§| ) - (5.2c ) could be written in the same manner as the expressions (5.12) - (5.14) showing potential WW- like 
expressions for the twist-2 distributions /^°n2 ^' x)- 

For the symmetric part of the Compton amplitude we perform the same calculational steps as in the anti- 

symmetric part, namely: 



• Separation of T'p/C"''— contractions; this also includes the trace terms. 

• Insertion oi V = P+ t + tt + tt Zk into the Compton amplitude (4.1C) 



• Definition of the structure functions F° „ using (5.5) 



Partial integration with respect to t by formula (5.9) 
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• Projection onto the coUinear part. 
After partial integration the symmetric part of the Compton ampHtude has the following form: 

r-l 



-2 



1 e + < - le 



1 



qP+ 



gfj.i> Qp /oo ~ {g^pQi^ + 9vp1pi) Fqo + {gpipP+v + 9ypP+p) tPoo 

+ {9p.pT^u + QupT^p) t /'Co + {gppT^v + g^pT^p.) t 
{q^P+, + q,P+,) (fS, - fSo) - P+pP+u {2tFSo -tfSo) 

+ {qp7T. + quTTp) (^F^o - /fo) - {P+P^u + P+uTTp) (t F^^ - t f^, 

p+. 



(5.21) 



+ {qpn. + q.^p) [f^, - f^,^ - {P+^n, + P+,n^) (t F^^ - t f^. 



2t i^20 t /2c 



(2tFS, 



- [tT,,^, + TT,^^) (2t F^^ ~ t /fi) 



^ap 



The collinear part of which contains only functions F^^ and /qq, can be written as 



^{pu} 




Looking at the tensor structure of ( ^.22| ), the functions 

F^(i;,7,x) - 2t./o"o(t;,7,x) , 
appear as natural structure functions and obey a Callan-Gross like relation 

F^(i;r7,x) = 2t.F?(i;,7,x) ■ 



P+„P+, 2t/°o(t;r;,x) (5.22) 

tFS^{t;r^,x)^]C^P . 

(5.23) 
(5.24) 

(5.25) 



Similar to Ref. § one observes that the remainder contributions in Eq. ( |5.22| ) are suppressed for large values of a 
property of the off-coUinear terms. One may see this, contracting the structures /C"'' with the respective tensors in 
front. As well known from forward scattering, the Callan-Gross relation receives corrections both from higher orders 
in the coupling constant and due to mass effects, see e.g. Ref. and therefore as well in the non- forward case. On 
the other hand, the Wandzura-Wilczek relation for geometric twist-2 turns out to be rigidly stable, see e.g. Ref. [P4| 



Including also 'Pp/C'^-contra ction s in the calculation will result in an additional tensor structure analogous to ( 5.17| ) 
with qp replaced by kp in Eq. (5.22). The corresponding structure functions will not be given in explicit form. 



VI. EVOLUTION EQUATIONS FOR THE DISTRIBUTION AMPLITUDES 

The scaling violations of the operator matrix elements and distribution amplitudes of the process considered are 
described by the renormalization group equations governing the ultra-violet behavior of the light-cone operators. 
The corresponding equations for the distribution amplitudes fa (z^, Z-, Zk) are c alled evolution equations, which we 
are going to discuss in x— and z— space. Since the flavor content of the operators (2/7) and (2^) has been suppressed 
in the preceding sections, we treat only the flavor non-singlet evolution equations as an example. The singlet evolution 
equations are of quite similar structure (see e.g. Ref. for earlier work see Rcfs. |33t). 
The non-singlet renormalization group equation for the twist-2 vector operator reads 



KiX, K2X; /i^) 



K [Ki, K2; K.-^^, K2) Uf^ [KiX, K2X; IJ, ) 



(6.1) 
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By contraction with it is obvious that the scalar twist-2 operator O*^^ „ ^0Qtw2 Qi^gyg exactly the same renor- 
malization group equation because multiplication with commutes with the differentiation on the left and with the 
integration on the right hand side. This gives the renormalization group equation for the scalar operator which on 
the light cone already is of twist-2: 



'-—rO(nix,K2x;iJ,^) = I DV 7 (/ti, K2; k'i, K2) O (k'iJ, Kji; /z^) . 
d^ JR2 



In the last equations the integration measure 

DV = d/i'i dK2 e (ki - k'i) e (k'i - K2) 9 (ki - 4) 6 (4 - K2) 



(6.2) 



(6.3) 



has been introduced. In Refs. ^ it is shown that the non-local anomalous dimension matrix 7 is invariant under 
translations and scale transformations, 



7 (ki, K2 ; k'ij '^2) = 7 - '«0, K2 - Ko ; 4 " ^0, 4 - Kq) 

= 7 (Aki, Ak2 ; Ak']^, AK2) 7 



(6.4) 



which reduces the number of independent variables of 7 by two. By first changing the variables from ki^2 to 
followed by a translation by k+ and a scaling by one derives the following form of the evolution kernel 7: 



7(ki,k2 ; 47«2) = 7(«:+,k; '«+,«^') 

= 7(0, k; 



(6.5) 



^7(0,1; 



K {wi,W2) 



with 



Wi 



and 



W2 



(6.6) 



The variables hi[ and Wi are connected by the following transformation 



1 - 



W2 



(6.7) 



It is therefore more natural to use wi and 'W2 as integration variables in the renormalization group equation (6.2) 
instead of k[ and The integration measures are related by 



(6.8) 



where and D\j include the suitable 6'-functions realizing the integration ranges of (6.3) and 



Wf. 



- dwi dw2 6* (1 + Wl — W2) {1 — Wl + W2) {1 + Wl + W2) {1 ~ Wl — W2) 



The measure TPw can be divided into two parts, because K (wi, W2) under the exchange Wi <-^- — wi, W2 
the following relations, cf. Ref. |^ : 

K {wi,W2) = K { — Wl,W2) = —K {wi,—W2) ■ 



Putting (|6.5|) and (3.S) into the renormalization group equation for the scalar operator results in 



/^^^j^O {kiX,K2X\^1^) 



TP'w K {wi,W2) O {^k'^x , K^x; fj?^ 



-W2 obeys 
(6.9) 

(6.10) 



For the explicit structure of K (wi , W2) see Refs. |, H, ||. The equation ( |6lo| ) wiU now be considered for the matrix 
elements of the scalar operator which are, according to equation ( 3.20| ), given by 



(p2, k\ O {kix, K2X, /i^) \pi) 



Ak-l. xP- v^a ( ~ 



JC^ {x,P2, k,pi) fa (i^xP^, KxP^, Kxk, fi^) . 



(6.11) 
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^From this one obtains directly the evolution equation for the distribution amplitudes fa in x-space: 
o d 



.12) 



Because we are interested in evolution equations in 2:-space, we perform a Fourier transformation of equation ( |6.12 ) 
The physically relevant transforms of fa are given by 



Jm Jr 2tt Ztt 

Carrying out these transformations one arrives at the following result: 



fa {KS:P+,KS:P^,Kik-^i^) .(6.13) 



fa{z+,z-^,zk,fi ) = , o ; o to 
d/i^ ./» zvr ./» 27r ./» 27r 



JR 27r Jjj 2tt Jjj, 27r 

X ( T)^Z' fa{z'_^,z'^,z'^;fl^) e-i^2>^S:{P+z'^+P^zL+kz',) ^iw^nxP^ 



/ D^z' /q (z^, z'_,z'k] H^) 6 {z+ - W2Z'^) S (z_ - W2z'_ + Wi) 6 {Zk - W2z'f,) 

JR3 



(6.14) 



with D'^z' given by (3.23) and D^z' defined by 

DV = 2dz^dz^ e'(l-z^ + z^)e'(l + z^- z^) 61(1- z^-z^) 61(1 + z^ + z^) 
We introduce the evolution kernel 



r [z+,z^,Zj^,z_) = - K [z_—-z^, — 

\z+\ \ z; z; 



which leads to the this evolution equation 



fa{z+,Z_,Zk 



= / dV r(z+,. 



Z— , Z I , Z_ ) fa ( Z I , Z_ , Zfc ; ^ 

z+ 



(6.15) 



(6.16) 



Let us point to the remarkable fact that the variable z^ connected to the meson momentum k only appears as a 
parameter in fa and is not contained in the evolution kernel F. The same observation has been made in Ref. |l^ 
recently in the case of diffractive scattering, where the parameters rj or xp behave in the same way. In so far some of 
the scaling variables of a problem, in the present case the variables Xij pla-y another role than others, as here ^ and 
r], which interfere with the evolution. 

This evolution equation is a fundamental equation because it describes the evolution of the triple- valued distribution 
amplitudes fa {z+, Z-, Zk) in z— sp ace. T hese amplitudes are the basic objects for the construction of the structure 
functions Fa,F^ and F*'' in (B.29) and (3.30). They are also used in the definition of the single-valued functions 



Fnin2 {tiViX) in The scaling violations of these functions are obtained solving Eq. (6.16) and inserting the 

functions fa into Eqs. (|3^ , ^^ . 

It is also possible to obtain another evolution equation for /qq (t; -q, x) in the variable which is compatible with the 
former equation. This single- variable evolution equation governs the evolution of the structure functions contained in 
the coUinear part of the Compton amplitude. 

To begin with, we first show that the distribution amplitude /qq (i; 77, x) given by 

(6.17) 



(6.18) 



/oo itViX) = j <^z^ j dzfc {t-T]Z- - xzk, z_, Zk) 
has another representation obtained as 

JC'^{p)fSo{t;V,x)^ J ^^^^e"^'''^^ {p2,k\0{~Ki,Ki) 
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The constraints 



xk = X iP+ 



(6.19) 
(6.20) 



appearing in the former equation are the scahng relations ( |l.7| ) in x— space. Using the representation (3.25) under 
these constraints leads to the result (6.17). 

To derive the single- variable evolution equation we form matrix elements of equation ( 6.10| ), 



/^~3~^(P2, ^1 O {-KX, KX) \pi) 



d/^2 



(6.21) 



D^W K{wi,W2) e^i"*-^- {P2,k\ O {-W2KX,W2Kx)\pi) 



xP- —f] xP-i^- ; xk—x xP-^ 



and perform the Fourier transformation in the variable kxP^ according to ( 6.1q ). The direct calculation leads to 

d 



with the evolution kernel 



-f{t,t';7]) = / dw2 



1 ~fw2t'-t 



,W2 



(6.22) 



(6.23) 



Like the evolution kernel T z_; z^, z^L) , also j{t,t';rj) does not depend on any fc— dependent variables like x or 
Zk being related to the meson momentum. 



VII. GENERALIZATION TO AN ARBITRARY NUMBER OF OUTGOING MESONS 



In this section we summarize the generic properties of the results obtained in the preceding sections and extend it 
to an arbitrary number of outgoing mesons. Generally, one may state that all the above results remain valid under 
slight modifications if two or more outgoing scalar mesons are present in the process, 

7r(gi) + H(pi) 7*(g2) + H(p2) + M(A:i) + ■ • • + M(fc„) , (7.1) 

as shown in Fig. |l|. To fix the kinematic domain of this process, all meson momenta ki are connected to different 




scaling variables Xi defined by 



FIG. 1: Process with n outgoing mesons 



where P-|_ and P_ are obviously given by 



Xr = 



(7.2) 



P± = ip2 + ^h) ±Pi , 



(7.3) 
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and ^ and rj are introduced as in (^]^). In order to compute the twist-2 part of the Compton amphtude 



RT 



J. 



(f)*(-|) 



S 



Pi, Si 



(7.4) 



for the general process (7.1) one apphes the same approximations to the operator T^i, as in Sections. ^ and [II. Thus 
one uses the approximation (2.12) and apphes the twist-2 projection (3.6). 

The technique used in Section III to construct the matrix elements (^p2,k \ O^*^-* \ pi) can be carried out for an 
arbitrary number n of scalar mesons where the additional meson momenta {k2, ■■■,kn} enlarge the set of kinematic 
factors lC°'{x,p). However, these factors are easy to guess and their number is given by 



r 2-+1 + forn > 1 

scalar ~ i o -C r\ 

2 lor n = 



(7.5) 



for the scalar matrix element. This formula also reproduces the number of kinematic factors in the ordinary non- 
forward case: the Dirac- and Pauli-structures. 

Having all kinematic factors /C" at hand, one introduces the related structure functions fa and writes down the 
following decomposition of the scalar matrix element 

{p2,ki,...,kk\0{-KX,Kx)\pi) = ^ ICa{x,p2,ki,...,kn,pi)fa{i^S:P+,KxP-,Kxki,...,Kxkn,PiPj;fi'^). 

a=l 

With this representation one goes through the same steps of the calculation as in the preceding sections, namely 

• Fourier transformation of fa 

• Application of the twist-2 projector 

• Computation of the Compton amplitude. 



The result is again of the form ( [4.5D with a larger z— space and V given by 

n 



(7.6) 



In this sense, the fo rm ( 4.5) of the Compton amplitude is a generic result holding for a large class of processes. The 
more explicit form (4.1C )Ts generalized by replacing X ^fe by ^ Xi ^ki ■ 

It is even possible to interpret the Wandzura-Wilczek and Callan-Gross relations obtained in Section ^ as generic 
properties of the coUinear parts of these processes. Making the substitution 



P -P+ t + TT Z_ + ^ TTi Zfc, with TTi = fcj - XiP+ 



(7.7) 



and projecting onto the coUincar part one finds again the relations 



/•sign(t) 

Fa (^;?7,Xi,---,Xn) = 2t- (i;?7,xi,...,X„). 



■ ,Xn) 



(7.8) 
(7.9) 



Looking at the derivation of the evolution equation in Section VI , it is not difficult to find the appropriate generalization 
to an arbitrary number of mesons (k = (fci, . . . , fc„)): 



d 



1+ 



(7.10) 



For n ~ Q, this general evolution equation is reproducing the fiavor non-singlet part of the evolution equation given in 
Rcf. Q for the ordinary non-forward scattering. As in the case of one meson in Section VI the variables z^^ connected 
to the meson momenta only appear as parameters and do not contribute to the evolution kernel. 



■ K \ z_— — z- , -r 
z'+ 



(7.11) 



20 



Only the number of mesons is relevant for the structure of this kernel. The single-variable evolution equation is of 
the same form as in the preceding section. One only has to enlarge the number of scaling parameters Xi- 

f^^^fSo{t;V,Xu---,Xn,fi^) = J ^dt' j{t,t';Tj) fSo{t';v,Xi,---,Xn,fi^) ■ (7.12) 



VIII. CONCLUSIONS 



We studied the structure of the virtual Compton amplitude for deep-inelastic non-forward scattering 7*(gi)+H(pi) — > 
7* (92) + H(p2) + M(fc) at the level of the twist-2 contributions in lowest order in QCD in the massless limit. In the 
extended Bjorken region, i.e., {{qP+),—q^} — > 00 with —q^/{qPj^),{qP^)/{qP+) and {kP^)/{qP+) kept fixed, the 
twist-2 contributions to the Compton amplitude were calculated using the non-local operator product expansion 
for general spin states. In this approximation the Compton amplitude consists of five kinematically independent 
parts which in the limit fc ^ reduce to the well known Dirac- and Pauli-type amplitudes. A decomposition of 
the Compton amplitude was performed with respect to the helicity states of both (virtual) photons. In complete 
analogy the (electromagnetic) gauge invariance of the non-local light-cone expansion holds at the level of the 5*- 
matrix since the fact that the leptonic currents are conserved. Due to this, only those contributions in the Compton 
amplitude are projected out, which obey gauge invariance. Integral relations generalizing the Callan-Gross and 
Wandzura-Wilczek relations for unpolarized and polarized forward-scattering are derived by reduction to the coUinear 
parts of the Compton amplitude and, thereby, reducing the triple-valued distribution amplitudes to one-valued ones 
(for {qP-)/{qP+) and {kP-)/{qP+) fixed). In this connection attention has been drawn to the difference between 
geometric and dynamic WW-relations being related to different notions of twist. The evolution kernels of these 
distribution amplitudes are obtained from the (well-known) non-local anomalous dimensions of the (scalar) twist-2 
light-ray operators 0*"^(Kii, K2i) and 0^*^^(Kii, K2i); they are independent of the meson momentum k. These 
results show that deeply virtual Compton scattering off nucleons in the case of additional meson production behaves 
quite similar to the case where mesons are absent. Both the basic structural relations as well as the scaling violations 
are the same in both cases. However, the structure of the Compton amplitude is different in general, however only 
with corrections of 0{1/^) or less. 
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APPENDIX A: PROJECTION OF THE TWIST-2 OPERATOR ONTO THE LIGHT-CONE 



This appendix is devoted to the derivation of the results ( 3.32| ) and ( 3.33 ) from the off-cone twist-2 non-local quark- 
antiquark operators obtained in Ref. |l8| . 

As has been shown there the nonlocal quark-antiquark operator of geometric twist-2 has the following structure 
((■07m V') (9) = Jd^ exp{i(ga;)}-(/;(x)7^?/;(-x)): 



d^g 
(2^ 



dr 



(3-Hg9g)(5^-irq^a;„ ^2(9^2;) 



(8.1) 



(3 + qdq) ((4 + qdq) itq^x^ - ^iir)^ {q^xf^x^ + x^q^qo.)) + l{itfq^q^Xa^x^] Hs{q\rx) 



Its n— th moment is given by 



1 



{h^.i^){q)\ 5^ hl{q\x) - q>^X^ hl_M^) 



{n + lf J (27r)4 

2x^'qa hl_^{q\x) - (x^x^g^ + qt^q^x^) hl_^{q\x) + | q^'xcX^q^ hl_^{q\x) 



(8.2) 
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Here, for notational simplicity we used the following abbreviations: 

1/2-1/ 



(8.3) 

The relation between the non-local and the local operators, Eqs. ( ^.1[ ) and ( ^.2D , off the light-cone is obtaine d by 
observing that the Bessel functions are generating functions of the Gegenbauer polynomials (see, e.g., Ref. [p4[ , 
Eq. II.5. 13.1.3): 



^ (2i^) 



(8.5) 



where (2v)n = 2i^(2i^ -f 1) . . . (2i^ + n - 1) = r(n + 'lv)lV{2v) is the Pochhammer symbol. 

The projection onto the light-cone is obtained most easily by first considering the local operators. Because of the 
series expansion of the Gegenbauer polynomials (see, e.g., Ref. Appendix 11.11), 



/£ = 



k\{n-'lk)\ ^ ' 



.6) 



one observes that from the expression (^^) on the light-cone, — 0, only the term with the highest power, i.e., for 
A; = 0, survives: 



{n + iy - 1)! 
nliu - 1)1 



■7) 



Using these results in the expression (3.2) we obtain: 



1 



n+1 
1 



{n + 1)2 J {27t) 



^ (V^TmV^) (9)£a [q^ n (gi)"-i -|- q^n{n- 1) (g^)"-^}. (8., 



Now, using 



1 



n -t- 1 



At r" and 



we are able to re-sum over n according to 



(n + l)2 
Ik)" 



dr r" In r 



n=0 



(8.9) 



(8.10) 



There are two options of doing this. In the first instance we may replac e in qg in Eq. (B.8) by the derivative 9^ acting 
on the exponential exp{iK(gi)}. Thi s way one retains the expression ( |3.5| ) of the non-local twist-2 operator on the 
light-cone from which the expr ession ( |3.32) has been derived. On the other hand, after taking matrix elements of (8.1) 
and observing t he d efinitions ( 3.29| ) and ( p. 30 ) of the distribution amplitudes Fa{z) and F*''(z), one obtains exactly 
the expression ( |3.32 ). Analogous results hold for the axial vector case (3.33). 

The same result could have been obtained also using the Poisson integral for the Bessel functions (cf., Ref. 
Eq. II.7.12.7), 



r(i.+ i)J,(z) = ^(^ 



d^(l-^2^I^-l/2 gitz 



for Re z/ > — ^, in order to express the functions ( |8.3D on the fight-cone by 

1 



H^{q\x) 



m Jo 



dA[A(l - A) 



''-I gi-^ (<?£) 



.11) 



.12) 



Then, after shifting the homogeneous derivations qdq in the expression (8.1) to the right and interpreting it as Xdx 
acting on the exponential, some partial integrations with respect to A can be performed which, finally, lead again to 
the expression (3.5) and (3.32), respectively. 
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APPENDIX B: HELICITY PROJECTIONS AND CURRENT CONSERVATION 



In this appendix we construct the hehcity projections of the Compton amphtude generahzing the results of Ref. 0| 
to the present case. We start with the construction of the hehcity basis of the two virtual photons 7^ and 72 . To 
simplify this construction, we choose the Breit frame, in which the relevant momenta read: 

(9.1) 

(9.2) 
(9.3) 
(9.4) 



P- = (0;0,0,P-3) , 
k = (fco; fci, fc2, fcs) 
q = (go; 91,0,53) , 



where /j, is introduced as a mass scale of the hadronic momentum P+ with = /i^. To define the helicity basis we 
introduce the two reference vectors 



no = (1; 0,0,0) 
n2 = (0; 0,1,0) 



The polarization vectors of the photons 7J' and 72 are then given by 



-1 



1 
1 



(9.5) 
(9.6) 



(9.7) 
(9.8) 



■'If, 



= 

=■2 _ 



-2 



"-2m 1 

f a /3 7 

^V21 

^ a 13 -r 

AT £^0/37 "-0 "2 I2 J 
^V22 

^ (g2M 92 -"-0 - ?^0A^ 92 -92) 



The normalization factors are given by 

.1/2 



iVoi = 



iVo2 = 



p2 



4i/ 



iV21 = - 



/;2 u2 P2 

l + ^(e-^)-^ 



Af2 



^31 = 



32 



1/3/2 

M 
i.3/2 

M 



;/2 /t2 P2 



V 


-rif-Pl 






-?/)H 


Ai^P^ 
16j/3 


V 


^^f-Pl 




2v^ 


h?/)^ 


^2p2 

16j/3 



(9.9) 
9.10) 

9.11) 

9.12) 
9.13) 

9.14) 
9.15) 
9.16) 
9.17) 
9.18) 
9.19) 
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which follows from the relations 





= V , 


(9.20) 




= rjv , 


(9.21) 


(qq) 


= -i^ , 


(9.22) 






(9.23) 



(see definitions (1.7) and note that ng = P+/ in the Breit frame) and conservation of momentum 



qi = q 



q2 = q 



2 



(9.24) 
(9.25) 



(9.26) 



The polarization vectors obey the following normalization condition, 

^a/i ~ "a '-'ab j 

with Sa = —1 for a = 0, 1, 2 and Sa = 1 ior a = 3. We now use this helicity basis to compute all matrix elements 

Jki^el''T';:^el'' with fc,Ze {0,1,2,3} (9.27) 
The result of the straightforward calculation is 

2 1 



00 



^ ^/\e^\ 



Tfi « 2 D^z - 



3 

1 q.r 



Q2 + ie (Q2+ig 



• \2 



' 22 



' 33 



2 / D^z I - 
2 



1 



q-V 



Q2+ie (Q2+ig)2 



FaqpK.^' 



Faiz+,z^,Zk) qpK."" 



(9.28) 
(9.29) 

(9.30) 
(9.31) 



-tF -tF -tF jl 
' OlJ ' 10' ' 02i ' 2 



03 



30 



' 12: ' 21 



-jF -jF -jF jF 
' 13' '31' ' 23' ' 32 



O 

2 



o 
o 



1 



1 + ^(^ + ^)1 



D3zF,(z+,Z_,Zfe)qp/C'^'' 



D3zFa(z+,Z_,Zfe)qp/C"^ 



1 



(9.32) 
(9.33) 
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for the symmetric part. The helicity projections of the trace terms are all of order l/v and therefore not given in 
explicit form. The same calculation is is carried out for the antisymmetric part: 



' 00 ' ' 11 — ^ 



' 22 I ' 33 — 



-rF5 -rF5 -rF5 -rFb _ 
' 01 I ' 10 ' ' 02 J ' 20 ^ ^ 



jFb jF5 _ 
' 03 ' ' 30 — 



-rFb -tF^ — n 

' 13 ' ' 31 ' ' 23 J ' 32 — 



1 

1 

1 

1 

7^ 



Tf2^ « 2 / D^z ( -1- - F^^ q^K^'^" (9.34) 



le 



(8' 



Here, we have only kept terms contributing to the highest power in because all other terms vanish in the limit 
V ^ Qo. Terms proportional to in the normalization factors have been neglected, because they do not contribute 
to the highest power in v. The amplitudes T^q,T|3 and T^3,T|q are a priori not of order \/v^ because gp/C"'' is of 
order v. But since the integrals 

/ D^z i^(z+,z_,Zfc) = (9.36) 

vanish, these amplitudes are also identical to zero. 

In the above only the contractions of the helicity vectors with the Compton amplitude were considered. For the 
physical process, however, the corresponding projections for the leptonic tensors L^^'^ have to be considered as well to 
see, which terms contribute to the physical S'-matrix. Due to the fact that 

L^li = Ll,<l2 ^ , (9.37) 

holds all remaining terms in the projections Tqa; and T^o are annihilated. 

In leading order in v only the amplitudes Tfj^, Tfj and Ti2, give a non- vanishing contribution in the extended 
Bjorken region, whereas other terms T^;"* for k,l — 1, 2, 3 are suppressed at least in 0{1/ The explicit calculation 
also shows that 

Tfi = T|; and Tf^^ = -T^i' (9.38) 

in leading order. Only two of the sixteen amplitudes are relevant for — > 00. Similar results have been obtained in 
Ref. [061. 
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